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The hyperfme structure (HFS) of a bound electron is modified by the self-interaction of the elec- 
tron with its own radiation field. This effect is known as the self-energy correction. In this work, 
we discuss the evaluation of higher-order self-energy corrections to the HFS of bound P states. 
These are expressed in a semi-analytic expansion involving powers of Za and \n(Za), where Z 
is the nuclear charge number and a is the fine-structure constant. We find that the correction 
of relative order a(Za) 2 involves only a single logarithm \n(Za) for P1/2 states [but no term of 
order a (Za) 2 \n 2 (Za)], whereas for P 3 / 2 states, even the single logarithm vanishes. By a Foldy- 
Wouthuysen transformation, we identify a nuclear-spin dependent correction to the electron's tran- 
sition current, which contributes to the HFS of P states. A comparison of the obtained analytic 
results to a numerical approach is made. 

PACS numbers: 12.20.Ds, 31.30.Jv, 31.15.-p, 06.20.Jr 
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I. INTRODUCTION 

Throughout the history of quantum electrodynamic 
(QED) calculations of atomic properties (notably, en- 
ergy shifts) , two approaches have mutually inspired each 
other, namely, the analytic and the numerical ansatz. 
The necessity for employing both methods is easily seen 
when one considers the range of coupling constants Za 
which are relevant for hydrogenlike ions (here, Z is the 
nuclear charge number, and a is the fine-structure con- 
stant). For low Z, the parameter Za is very small, and 
an expansion of the fermionic propagators in powers of 
Za appears reasonable. For high Z , the parameter Za 
approaches unity, and the mentioned expansion is not 
practical [l[ . Consequently, starting with the accurate in- 
vestigations on the spectrum of highly charged ions @, Q , 
there has been tremendous activity over the past few 
decades regarding the accurate description of energy lev- 
els of heavy ions. These are complemented by technically 
demanding experiments Q. 

Some of the most important effects to be considered in 
QED bound-state calculations are so-called self-energy 
corrections, where a bound electron spontaneously emits 
and reabsorbs a virtual photon. Between the photon 
emission and absorption, other interactions — with the 
binding Coulomb field and, possibly, with other external 
fields — may occur. Even for low-Z ions, a direct expan- 
sion of the electron propagators in powers of Za is not 
universally possible: namely, the energy of the virtual 
photon must be large enough (larger than the scale of 
the atomic binding), so that the electron can essentially 
be regarded as a free particle in between the emission 
and absorption (perturbed only by a finite number of in- 



teractions with the Coulomb field). However, when the 
photon energy is small (commensurate with the atomic 
binding energy), this expansion is no longer possible. In 
this case, the expansion in powers of Za is achieved "im- 
plicitly," by observing that since the photon energy is so 
small, one may expand the currents at the photon emis- 
sion and absorption vertices in the long- wavelength limit, 
i.e. in terms of dipole interactions, quadrupole interac- 
tions, spin-dependent interactions etc. (see Ref. 0, 0]). 
The Foldy-Wouthuysen transformation 0] can be used 
in order to achieve a clear separation of the Hamilto- 
nian into a leading nonrclativistic term and relativistic 
corrections, and a decoupling of upper and lower com- 
ponents of the Dirac wave function to a specified order 
in Za is achieved. This transformation, together with 
a careful matching procedure needed in order to "join" 
the high-and low-energy parts, then leads to the analytic 
results traditionally used in order to describe the Lamb 
shift [a]8| and other effects such as the bound-electron g 
factor Q. 

The intricacies described above are responsible for the 
Za expansion not being a simple power expansion (Tay- 
lor series) in Za. The matching of high-and low-energy 
contributions at an intermediate photon energy scale 
commensurate with an overlapping parameter e leads to 
the appearance of logarithms of Zajsee also the illus- 
trating example in Appendix A of jTy] ) ■ The expansion 
thus is semi-analytic. For the Lamb shift in hydrogcnlikc 
systems, many terms have been calculated in the semi- 
analytic expansion [||, but the numerical approach was 
faced with tremendous problems, and no direct result had 
been calculated for hydrogen (Z = 1) up to the year 1999. 
At low Z , the experimental accuracy is orders of magni- 
tude higher than at high Z , and the comparison and the 
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implied mutual consistency check of the numerical and 
analytic calculations is most meaningful. Therefore, a 
calculation was carried out [Hj] which confirmed the con- 
sistency of both approaches and determined a nonpertur- 
bative remainder term which is beyond the sum of the 
known terms in the Za expansion. The nonperturbative 
remainder term amounts to roughly 28 kHz for the hydro- 
gen ground state Lamb shift, which is a numerically large 
effect as compared to the current experimental accuracy 
of about 33 Hz for the 15-25 transition [l^, EH- Similar 
calculations were done, and agreement of the a naly tic and 
numerical calculations was found, for P states [lif . Other 
effects studied both within the Za expansion and within 
the numerical approach include the bound-electron g fac- 
tor of a bound S state 0, [Hj]. Also, the self-energy cor- 
rection to the hyperfinc structure (HFS) was extensively 
studied for the S states during last decades, both within 
the numerical all-order approach [l6l . [TtI [l8j and within 
the Za expansion [13 , l2fl] . 

For the hyperfinc splitting of P states, however, in- 
vestigations of the self-energy corrections are much more 
scarce, both within the numerical as well as within the 
analytic approach. There are some quite recent all-order 
numerical calculations for the 2Pi/ 2 state [2l[ as well as 
for the 2P 3 / 2 level [13] (see also the latest paper [23j). 
but there are only few analytic results to compare to. If 
we denote the nonrclativistic Fermi energy by Ep, then 
the only known self-energy correction [24[ is that of order 
^Ep. This correction amounts to j^Ep for nPi/ 2 states 
and —-^Ep for nP 3 / 2 states. Because this correction is 
entirely due to the electron magnetic moment, we can 
refer to it as an anomalous magnetic moment correction 
to the HFS. 

In this work, in the order to address the current, some- 
what unsatisfactory status of theory, we calculate the 
self-energy correction to the HFS of P states up to or- 
der a(Za) 2 EF . This correction goes beyond the anoma- 
lous magnetic moment correction and is due to numer- 
ous effects. Indeed, the calculation constitutes quite a 
complicated problem, mainly for three reasons. First, 
the problem considered is a radiative correction under 
the influence of an additional external field, i.e., we have 
three fields to consider (the photon field, the Coulomb 
field and the nuclear magnetic field). Second, the angu- 
lar algebra is much more complicated than for reference 
S states [25j. The third difficulty is that for P states, the 
nonrclativistic limit of the hyperfinc interaction has to be 
treated much more carefully than for S states. Namely, 
we may anticipate that there exists a correction to the 
electron's transition current, caused by the hyperfine in- 
teraction, which contributes to the self-energy correction 
for P states, but vanishes after angular integration for 
S states. This correction to the current can be under- 
stood easily if one considers the coupling of the physical 
momentum p — e Ahf s of the electron to the vector po- 
tential of the quantized electromagnetic field [here, 
is the vector potential corresponding to the hyperfine in- 
teraction defined below in Eq. J2J)]. The correction to the 



current has not been considered in the previous investi- 
gations which dealt with S states [Til Hfl . 

We organize our investigation as follows. In Sec. [ill 
we present some known formulas needed for the descrip- 
tion of the HFS. In Sec. IIII1 we consider the Foldy- 
Wouthuysen transformation of the Hamiltonian, and of 
the transition current, and identify all terms relevant for 
the current investigation. In Sec. IIVI the logarithmic 
terms of order a(Za) 2 \n(Za)Ep are given special atten- 
tion. Their value is derived within a straightforward and 
concise analytic approach. We then continue, in Sec. [Vl 
to investigate the contribution of high-energy photons to 
the HFS of P states. Six effective operators are derived 
which are evaluated for general principal quantum num- 
ber of the reference state. The low-energy part is treated 
next (see Sec. EU). The vacuum-polarization correction 
is obtained in Sec. IVIII The results arc summarized in 
Sec. I Villi and conclusions are drawn in Sec. IIXI Natural 
units (h = c = eo = 1) are used throughout this paper. 



II. GENERAL FORMULAS 

We work in the nonrecoil limit of an infinitely heavy 
nucleus and ignore the mixing of 2Pi/ 2 and 2P 3 / 2 states 
due to the hyperfine splitting (this mixing is otherwise 
described in Sec. Ill C of Rcf. 26]). Under these as- 
sumptions, the relativistic magnetic dipolc interaction of 
the nuclear magnetic moment and an electron in is given 
by the Hamiltonian 

Hms = -ea- Ahfs(r) = |e| a ■ A ills (r) , (1) 

where the vector potential reads 

(2) 

•±71 I ~ 

so that 



a rsf\ 1 ^ x ? 



|e| a-(flxr) |e| g • (f x a) 
47r r J 47T r a 

Here, jl denotes the operator of the nuclear magnetic 
moment. In this paper, we will use the convention of 
labelling the relativistic operators by indices with lower- 
case letters and nonrclativistic Hamilton operators by in- 
dices with upper-case symbols. For future reference, we 
give the magnetic field corresponding to the vector po- 
tential ©, 



47r r 3 



(4) 



The operator -ffhfs acts in the space of the coupled 
electron-nucleus states 

\FMpIj} = J2C™£ l \IM)\jri, (5) 

where / and M is the nuclear spin and its projection, 
j and fi is the total electron angular momentum and its 
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projections, and F and Mp is the total momentum of the 
system and its projection (we denote the momentum pro- 
jection by fi instead of m in order to differentiate it from 
the electron mass m). With the help of the Wigner- 
Eckhart theorem, the expectation value of -ffhfs on the 
coupled wave functions can be reduced to a matrix ele- 
ment evaluated on the electronic wave functions only, 



(FMpIj | Hms I FMpIj) 

(2&Cj)) (FM F Ij 



(6) 
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l -^ l fUj) [F(F + l)-I{I+l)-j(j + l)} 
a9 ~f^ Uj) [F{F 1(1 jij + 1)] 



We have used |e| 2 = Ana and = giql\e\/ '(2m p ), where 
qn is the nuclear g factor and |e|/(2?7ip) is the nuclear 
magneton. The quantity £ e (j) depends on the electronic 
state only. Let denote the electronic state with total 
angular momentum j and angular momentum projection 
\i (we here suppress the orbital angular momentum in our 
notation for the electronic state). Then, with the index 
q denoting the vector component in the spherical basis, 
we have 
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(7) 



As evident from the second line in the above equation, 
£ e (j) does not depend on the actual value of the mo- 
mentum projection /i of the total angular momentum of 
the electron, as its dependence cancels between the nu- 
merator and denomcnator. In a similar way, the nuclear 
variables can be factorized out in evaluations of various 
corrections to the HFS, reducing the problem in hand to 
the evaluation of an expectation value of an operator on 
electronic state with a definite angular momentum pro- 
jection. In practical calculations, we always assume the 
angular momentum projection of the reference electron 
state to be = i, as in the third line of Eq. ([7]). 

In the nonrclativistic limit, the magnetic dipolc inter- 



action describing the HFS consists of three terms, 



#hfs =^—ft-h= ^— (h s + h D + h L ), (8a) 



h D = 



3m 2 

3 f{a ■ f) —a 
2 m 2 r 3 



(8b) 
(8c) 
(8d) 



Analogously to the relativistic case, the nuclear degrees 
of freedom in the expectation value of the operator Hrfs 
are factorized out, and the problem is reduced to an eval- 
uation of the matrix element of the purely electronic op- 
erator 



4 <T 6*(r) + 3{ *- f)fo - ao 



3m 2 



27rm 2 r 3 



(9) 



where f is the unity vector f/r. The nonrelativistic limit 
of £ e (j) is 



e e NR (j) = m^m 



(Za) 3 m 



\it\ n 3 (2 K + l)( K 2 -±) 



(10) 



where k = (— iy~ l+1 / 2 (2j + l) is the Dirac angular quan- 
tum number. We have defined £ so that it has dimension 
of mass (energy) and so that its normalization reproduces 
the characteristic 4/3 prefactor for the Fermi splitting of 
S states. For nPi/ 2 and nP 3 / 2 states, we have, respec- 
tively, 



e e NR (l)-e R (2P V2 )- 4( ^ )3m 



9 n 3 



4 NR (i) = er(2P 3/2 ) = ^ 



4 (Zafm 



45 n J 



(11a) 
(lib) 



Various corrections to the HFS can be conveniently 
expressed in terms of multiplicative corrections to the 
quantity £ e NR (j), 



(12) 



The corresponding corrections to the position of the HFS 
sublcvcls will then be 



8E HFS =E F 8S e (j) 



(13) 



where Ep (the Fermi energy) is the nonrelativistic limit 
of Eq. ([6]). In order to keep our notations concise, we 
define the normalization factor 



1 



1 



(14) 



which will be extensively used throughout the paper. An- 
other reason for our choice of the normalization of ho is 
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that we can use this operator as a perturbation Hamilto- 
nian for the ordinary nonrelativistic self-energy (with the 
correct physical dimension of mass/energy), in order to 
evaluate the relative correction to the Fermi energy, pro- 
vided we use a reference state with angular momentum 
projection /j = 1/2. 



III. FOLDY-WOUTHUYSEN 
TRANSFORMATION 

The Foldy-Wouthuysen transformation [3] is a conve- 
nient tool for obtaining the nonrelativistic expansion of 
the Dirac Hamiltonian in external fields. The idea is to 
construct such a unitary transformation of the original 
Hamiltonian that the transformed Hamiltonian does not 
couple the upper and the lower components of the Dirac 
wave function to a specified order in Za. In our case, we 
choose the starting Hamiltonian H t to be the sum of the 
Dirac Hamiltonian H re \, 



- - n % a 

H IC \ = a ■ p + p m . 

r 



and the relativistic HFS interaction operator, 



(15) 



(16) 



For the purpose of the present investigation we construct 
the Foldy-Wouthuysen transformation U that decouples 
the upper and the lower components of the Dirac wave 
function up to order (Za) 4 for contributions to the en- 
ergy and up to order (Za) 3 for contributions proportinal 
to the magnetic moment. Because the general paradigm 
of the Foldy-Wouthuysen transformation has been ex- 
tensively discussed in the literature for a large class of 
potentials [l?], HE] > we skip details of the derivation and 
just indicate the results. However, and this is an impor- 
tant point of the current paper, we must keep in mind 
that the Foldy-Wouthuysen transformation reads 



C/ = exp(i5), S = -i/?Odd(# t ) . 



(17) 



where Odd(if t ) represents the matrix of the odd compo- 
nents of H t in 4 x 4 spinor space, when the 4x4 ma- 
trix is broken up in 2 x 2 sub- matrices @, HI] . Because 
H t contains H^ s , the Foldy-Wouthuysen transformation 
chan ges as compared to an ordinary Lamb-shift calcula- 
tion [30fl . The transformed Hamiltonian H[ is 



Hi = UHtXJ- 1 = H 



FW 



H 



HFS 



(18a) 



where the first part i?FW does not depend on the nuclear 
moment. The second part is just the nonrelativistic HFS 
operator © , which we reproduce here on the basis of the 
Foldy-Wouthuysen transformation of the total relativis- 
tic Hamiltonian H t . The operator H-pw is a 4 x 4 matrix 
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\ 2m I r 
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nZa 



8m 3 2m 2 



6 3 (r) 



Am 2 7 



-E • L. 



(18b) 



Here, the E = 



are the 4x4 generalizations of 



a 
a 

the 2x2 Pauli matrices a. For the upper components 
of the wave function, JTfw can be replaced by the 2x2 
matrix 



Za 

-Hfw — * m + 7; 



V_ 

2m 



nZa 



8m 3 2m 2 



Za 
4 m 2 r 3 



a ■ L , 
(18c) 

where we identify the second and third term as the non- 
relativistic Schrodingcr Hamiltonian Anr- By contrast, 
for the lower components, the applicable 2x2 matrix is 
[up to order (Za) 2 ] 



Hpw 



(18d) 



p* Za 
2m r 

However, the transformation of the Hamiltonian is not 
the only effect of the Foldy-Wouthuysen transformation. 
In order to see that we also have to transform the transi- 
tion current of the electron, we remember that the char- 
acteristic integrand of a self-energy calculation is @, @] 

a 1 exp(ifc • r)— — a 1 exp(— ifc • r) 



H t - Et 



■ to 



V-' 



(19) 

where (to, k) is the four-momentum of the virtual photon, 
and E t is the total energy (corresponding to the Hamil- 
tonian H t ) of the relativistic reference state ip. In order 
to achieve a nonrelativistic expansion, we perform the 
Foldy-Wouthuysen transformation and write M. as 



M 



j' 



U(H t -Et + u)U 



-(ry 



(20) 



Here, \Utjj) is the nonrelativistic cigenkct correspond- 
ing to the Schrodiger-Pauli eigenstate (plus relativistic 
corrections and HFS-induced wave- function corrections), 
and we see how the transformed Hamiltonian H[ in the 
denominator is obtained [comparing to Eq. Q18a[) ]. The 
transformed current J\ which reads 



J 1 = Uacxp(ik ■ r)U 1 = c?exp(ifc 



P 



(21) 



contains the dipole current p/m and higher-order terms. 
These higher-order terms, in the absence of the HFS in- 
teraction, are listed in Ref. [3(|. In the presence of the 
HFS interaction, we find, however, an additional con- 
tribution to the current, due to the replacement p — » 
p — e for the momentum of the electron in the pres- 
ence of the hfs vector potential, 
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(22) 
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where fi is the unit vector p/\fl\. The contribution in- 
duced by Sj vanishes if radiative corrections are evalu- 
ated for the reference S states, but not for P states which 
are investigated here. In Eq. (f2"2")h we define the prefactor 
multiplying Sj in a way to be consistent with the normal- 
ization of the operator h in Eq. {8}. As discussed in the 
previous section, the nuclear degrees of freedom can be 
effectively separated out by assuming that the magnetic 
moment of the nucleus is pointing into the z direction. In 
this case, the correction to the current takes the following 
form in Cartesian coordinates, 



5 j -* Sj = - 



r-3 ^ ' 



(23) 



We note that, stricktly speaking, both the operator p/m 
on the right-hand side of Eq. ([2Tj) as well as the opera- 
tor on the right-hand side of Eq. (f2"2")) should carry a (3- 
matrix (30j . However, it can be replaced by unity when 
the current is applied to the upper components of the 
wave function, which are the only nonvanishing ones in 
the nonrelativistic approximation. 
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(28) 



The parameter e cancels at the end of the calculation, 
when the high-energy part is added. For the determina- 
tion of the logarithmic contribution it sufficient just to 
replace e by the electron mass m. To identify the correc- 
tion of first order in Hy^s, we expand the reference-state 
wave function as 



\<h) = \<l>) 



i 



i/HFS |<», (29) 



where the prime denotes the reduced Green function. 
The first-order pcrturbative correction of SE\ og is 



IV. LOGARITHMIC TERM 



In this section, we present a concise derivation of the 
logarithmic part of the self-energy contribution of order 
a 2 {Za) 2 Ep to the HFS of P states. To this end, we con- 
sider the perturbation of the nonrelativistic self-energy of 
the bound electron by the nonrelativistic hyperfine inter- 
action -ffHFS- Let -£/nr be the Schrodingcr Hamiltonian 
and Ht = -Hnr + -^hfs denote the total nonrelativistic 
Hamiltonian, whose the reference-state eigenvalue will be 
denoted by Et and the corresponding eigenfunction, by 
\4t)- The nonrelativistic self-energy correction of the 
state cf>T is 



8E=*L 
3tt 



dbJ i 



1 



m Et — {Ht + w) vn 



where e is a non-covariant frequency cutoff for the virtual 
photon. We use the expansion 



Et — (Ht 



Et — Hj 



O 



(25) 

which is valid in the domain u> G ((Za) 2 m,f) relevant 
to the calculation of the logarithm (note that in natural 
units, e has dimension of energy, or, equivalently, mass). 
The logarithmic part of the correction is generated by 
the integral 



SE } 



loir 



■ln[(Za)- 2 ] \{^\[p\ [Fhfs.P*]] 
f 



[p\ [i?NR,p 1 ]] 



EnV! — H] 



NR 



Hhfs 



(30) 



For P states, the second term in brackets vanishes be- 
cause \p l , [-Hnr, P 1 ]] is proportional to a Dirac S func- 
tion. Therefore, 



aAf 
3irm 2 

aAf 
3irm 2 



H(Za)- 2 ] {j\\[p\ [h,,p l ]]\j\) 



ln[(Zc 



V 2 /i 



(31) 



Here, \jh) is the Schrodinger-Pauli eigenstate \nPj) with 
angular momentum projection ^. For P states, we obtain 



(24) 
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n 2 — 1 (Zafm 




(32) 



The Kroncckcr symbol in the above equation implies that 
the matrix element vanishes for P3/2 states. The final 
result for the logarithmic part of the correction is 



1 a 



log 



(Za) 2 \n[(ZaY 



J '2 



(33) 



So, the self-energy correction to the HFS of P states 
can be conveniently expressed as 



a 

7T 



a o + {Za) 2 {a 2 i ln[(Za) 2 ]+a 20 } + 



(34) 
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where . . . denote the higher-order terms. As usual, the 
first index of counts the power of Za, and the second 
one indicates the power of the logarithm. 

The result ([33]) confirms the estimates of the logarith- 
mic coefficient 021 derived for the 2Pj states in Ref. [23| 
on the basis of an analysis of numerical data. Specifically, 
the values of a2i(2Pi/2) = —1.5 and 021(2/3/2) = 0.0 
were reported in that work, in full agreement with the 
analytical result of Eq. 



V. HIGH-ENERGY PART 

In this section we derive the part of the self-energy 
correction to order a(Za) 2 Ep induced by virtual pho- 
tons of high frequency, which is referred to as the high- 
energy part. It can be obtained from the Dirac-Coulomb 
Hamiltonian modified by the presence of the free-electron 
form factors Fi and F2 (for a derivation see, e.g., Chap. 7 

of mi), 



ffrad = c£-[p-eF 1 (V 2 ) A\ + /3 m + F 1 (V 2 ) V 

+ ^(V 2 )-^- (i^-E-pt-B 
2 m V 



(35) 



where V = —Za/r is the Coulomb potential. The form 
factors present in this Hamiltonian lead to various ra- 
diative corrections to the HFS, when A and B are re- 
placed by the vector potential and the magnetic field cor- 
responding to the hypcrfine interaction, respectively For 
S states, this procedure is described in detail in Ref. [25| . 

We find that the Hamiltonian (|3"5"|) induces six contri- 
butions of order a(Za) 2 Ep, 



(36) 



each of which will be addressed in turn in the following. 

The first correction C\ is induced by a term with -^2(0) 
in Eq. (|35p . namely 



F 2 (0)-^/?g.5=f- 
2m 2ir 



\e\m 

4-7T 



(1 ■ (h s + h d ) 



(37) 



Here f3 is the Dirac 7 matrix in the Dirac representation, 
^2(0) = ^r, and the vectors h s and hd are the 4x4 
generalizations of hs and ho, respectively, 



hd = 



3f(S-f) -S 



(38a) 
(38b) 



2 m 2 r 3 

The corresponding correction C\ is 

C * = l^ Qh l^(^,o + hd fi )\j\) R , (39) 



where h s $ and hd.o are the z components of the Hamil- 
tonian operators defined in Eq. (|38|) . We note that C\ 
contains the leading form-factor contribution of order a. 
To derive the next-order a(Za) 2 correction, one has to 
evaluate the matrix element with the relativistic (Dirac) 
wave functions (which have to be expanded in powers 
of Za beforehand in order to escape divergences due to 
higher-order terms) . By the index R, we denote the ma- 
trix elements evaluated on the relativistic wave functions. 

The second correction (C2) is an F' 2 correction to the 
effective potential ([37jl . i.e., 



F'M ^ /9 V 2 3 B 



a 

12tT 



4tt 



■ Pit- {v 2 (h s + h d )} 



(40) 



where we have used F^O) — jf-. Up to the order 
a (Za) 2 Ef, we can approximate the relativistic oper- 
ators h s and hd by their nonrelativistic counterparts hs 
and hp, replace the (3 matrix by unity, and write 



r oAf / • 1 
C2 = 12^ V 2 



V 2 (h s ,o + h D<0 ) 



(41) 



to be evaluated on the nonrelativistic wave functions, 

The third correction C3 is given by the term with F[ 
in Eq. ([35]) . namely 



V^hfi 
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where e is a noncovariant low-energy photon cut-off in 
the slope of the form factor F\. We can formulate C3 
nonrelativistically, 



C, = 



aN 
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(42) 



The forth contribution to the high-energy part is 
a second-order perturbative correction induced by the 
form-factor correction to the Coulomb potential V, 



: V 4 = Fi(V)-l 



1 


V 






/TON 


11" 


In — 






\2eJ 


+ 24 



V 2 

— v + 

to 



The corresponding correction is 

v,- 1 




-EWr — #NR 



(43) 



(44) 



where we have done a Foldy-Wouthuysen transformation 
on the propagator denominator in order to obtain the 
Schrodingcr Hamiltonian 7?nr and ignored higher-order 
(in Za) terms. Because V4 is proportional to a Dirac <5 
function, the correction C4 vanishes for P states. 
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The next contribution is the second-order pcrturbative 
correction induced by the relativistic hypcrfine potential 
-Hhfs as given in Eq. fl} and the following term in Eq. |35|) 

-V 56 = F 2 (0)—i^-E = -i- 7'Vy, (45) 

7r Zm 4irm 

where E is the electric field generated by the Coulomb 
potential V. The total correction then is 



Hi 



his 



Iff, 



his 



(46) 



It is conveniently splitted into two parts, C5 and Cg, as 
will be discussed below. We note that the relativistic 
HFS interaction 7J hfs couples the upper and the lower 
components of the wave function, as does V56, so that 
the second-order matrix element has to be evaluated on 
the relativistic wave functions (which is indicated by the 
index R). 

Let us consider the Foldy-Wouthuysen transformation 
of the numerator of the expression (|46[) very carefully. 
We write this numerator, employing a transformation hi, 
as 



hiV 56 hi~ 



1 



U{E Te \ — H re \)U 



UH Ms W 



(47) 

where hi, in contrast to U, is the Foldy-Wouthuyscn 
transformation that diagonalizcs the plain Dirac Hamil- 
tonian (fT5"|) without the HFS interaction (the state 
\<f>) after the transformation is just the nonrelativis- 
tic Schrodinger-Pauli eigenstate). In particular hi is 
obtained from Eq. (fl"7|) by the replacement H t — > 
H xe \. It has been shown in Ref. [3fJ that the Foldy- 
Wouthuysen transformation, when applied to a "third- 
party" operator — such as the electron transition current 
operator a exp(ifc-r) — leaves the leading term intact, and 
gives rise to higher-order corrections. We obtain in the 
case of V56 , 



1 



a 



8to 2 



v*v + — 



47r m 2 r 



Za 



a-L 



(48) 



V 2 ^, multiplied by 



Note that when we add the term ^-5 
the prefactor a/n, to the potential V4, then we obtain 
the effective one-loop Lamb shift potential, 



1 



8m 2 



;V V = 



a 
3^ 



(Za) 





5" 








+ 6 



m 



(49) 

which is useful when the entire formalism is applied to S 
states [see Eq. (21) of Ref. [Ill]- Finally, and somewhat 
surprisingly, the transformation hi ^ U applied to the 
relativistic hyperfine interaction i?hfs leads to 



UHhisU 



Hi 



hfs 



HFS 



(50) 



We retain the original relativistic HFS potential as 
the leading term and obtain the full nonrelativistic in- 
teraction i?HFS as an additional term, as well as higher- 
order terms which we can ignore. The Hamiltonian in the 
propagator denominator is transformed with the help of 

hiH^hi- 1 = H FW , (51) 

where i/FW is the 4x4 Foldy-Wouthuyscn Hamiltonian 
given in Eq. (|18b[) . which breaks up into the upper and 
lower components given in Eqs. (pel) and ipd"]) . respec- 
tively. From the upper components and the third term 
on the right-hand side of (|48p. we have 



a-L 



HFS 



\H 



HFS I 




Za 



1 



(52) 

The correction Cq is obtained when lower components of 
the Green function are selected. The bound-state energy 
can be approximated by m, and the energy of the virtual 
state as —m [see Eq. (fT8d| ], and we obtain from the first 
term on the right-hand side of 



C R = — 



2a ((j)\V 5 6 ^ H h fs I 4>) 



7T 

aN 

7T 



1 |Hhfs| 



(53) 



24 



^■,1/2 



GO 



4 



Za 



(54) 

One might wonder what would have happened if we had 
done the Foldy-Wouthuysen transformation in Eq. (|47j) 
with the full Foldy-Wouthuysen operator U, which also 
diagonalizcs the HFS interaction. In that case, we would 
not have obtained the term Hms on the right-hand side 
of (j5fj)) , but the transformation of the potential V56 would 
have yielded an additional term, proportional to the nu- 
clear magnetic moment. 

In order to check the above derivation of the sum of the 
C5 and Cq corrections, we evaluate Eq. |46|) in a differ- 
ent way, by using generalized virial relations (32l . [33[ for 
the Dirac equation. First, we introduce the relativistic 
perturbed wave function 5(j> as 



1 



[r x a 



(55) 



We note that the operator ' rxa j° i n the above expres- 
sion is the electronic part of -ffhfs after the separation of 
the nuclear degrees of freedom, sec Eq. ([7]). Performing 
the angular integration in Eq. (|46p . we obtain (for an 
arbitrary reference state) 

7a 2 KT r°° 

C 5 +C 6 =—^ r J o dr[g(r)Sf(r)+f(r)6g(r)}, 

(56) 
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where g and / are the upper and the lower components 
of the (relativistic) reference-state wave function, respec- 
tively, and Sg and Sf arc those of the (diagonal in k part 
of the) relativistic perturbed wave function 8<p. The per- 
turbed wave-function components 5g and Sf are known 
in closed analytical form [32I [33| ■ The radial integral in 
Eq. (|56p diverges for plain Dirac wave functions because 
of logarithmic singularities induced by higher-order terms 
in the Za expansion of the integrand. One first has to 
expand the wave functions in Za and then perform the 
integration. The result obtained in this way coincides 
with the one derived from Eqs. (|52"|) and J5T 



We now summarize the high-energy corrections, 



and those for nP$/2 states, 



ci(f) 



4+(*«) a (— - 3 



a !7 \2^Pl-}) 



c 3 (f) = o, 



C 5 (l) = >) 2 f-l-i 



360 An 10n 2 



7 



960 32 n 20n' 



(59a) 

(59b) 
(59c) 

(59d) 



A'' 



V 2 (ftsyj + 



3 2 1 1 






11" 












V 2 h 







Za 



a-L 



(57a) 
(57b) 
(57c) 




(57d) 



C 6 (|) = -(Za) 2 

TT 



3n 2 - 2 
10 n 2 



(59e) 



The net result for the high-energy part of the self-energy 
correction to the HFS of P states is 



TT A 



-2- 



19 7 5 
144 + 8^ + 16n 2 



1 / m 
— la I — - 

V2e 



(60) 



7T 



1 

24 



1 

60 





Za 




m 3 r 4 



(57e) 



^ H (l) = -' 



-- + (Zc) 2 



109 



2880 32 n 30n 2 



11 1 

+ 



(61) 



An evaluation leads to the following results for nP\/2 
states 



VI. LOW-ENERGY PART 



1 , 2 /13 3 7 

- + (Za) 2 1 - 

4 v ' \A8 8n 16n 2 



C-2 (3) =-(Za) — — 
z 7T 3n z 



C 3 (i) =-(Za) 2 ' 



2 In 



m\ 11 



2eJ 12 



C 6 (|)=-"(Za) 



a,„ , 2 3n 2 - 2 
20 n 2 ' 



227 1 3 

In 5 n 2 / ' 



(58a) 
(58b) 
(58c) 
(58d) 
(58e) 



In this section we derive the part of the self-energy cor- 
rection to order a(Za) 2 Ep induced by virtual photons of 
low frequency, which is referred to as the low-energy part 
in the following. In order to keep the notation concise, 
we suppress the indication of the reference state in 
the subsequent formulas; it is assumed that all matrix 
elements in the low-energy part are evaluated with this 
reference state. 

The low-energy contribution is conveniently separated 
into four parts that can be interpreted as corrections to 
the Hamiltonian, to the reference-state wave function, to 
the reference-state energy, and to the current. The cor- 
rection to the Hamiltonian is expanded first in Za, then 
in the overlapping parameter e. The standard procedure 
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(see, e.g., Refs. yields 



_L , .. 2aN , 



3tt 



P' 



1 



-h - 



1 



m Efiji — Pnr — w Pnr — #nr — urn 



2aM 
3irm 2 



In 



(2a) 2 m 



([p i ) N,P i ]])+(p 2 /io) 



+ -(Z«) 2 %(j). 



(62) 



Here, PhU) is a Bethe-logarithm type correction which 
needs to be evaluated numerically. The correction to the 
wave function is expanded into logarithmic and nonlog- 
arithmic parts as follows, 



1 



\ 171 PnR — ^NR — LU m \ -Bnr — #nr 
2aM 



3nm 2 



In 



(Za) 2 m 



(f) (h ) - (fh ) 



2([p\[(H mi -E NR ),p i }} 



+ ft), 

7T 



1 



E 



NR 



H 



Mi 



h 



(63) 



where we have used commutation relations. The correc- 
tion to the energy is 



37r J \m\ Enr - Hnr -oj J m 
o x ' 



2a 

37T771 2 



hi 



(Za) 2 m 



(p 2 ) + -{Zaf (3 E (J) ■ (64) 



The low-energy correction due to the nuclear-spin depen- 
dent current is 



^j(j) = -r— I dLuco 1 



3tt 



1 



m E — H — uj m 



sum of all discussed corrections gives the low-energy part 



37TTO 2 



hi 



(Za) 2 m 



{\p\ [ho, P 1 ]]) 



+ 2([p\\H NK: p l 



1 



(E-H 



-2{Za 2 



1 , / m \ , 



■(ZafpU), (66) 



where 



m = Mi) + /Mi) + Mi) + #/(i) • (67) 



In Eq. (f6"6"| . we took into account that 
(jp 1 , [H-NR,p 1 ]] ^e~h)' vanishes and the matrix 

element ([p 1 , [ho, p 1 )]) is evaluated in Eq. (132)) . It can 
be immediately seen that the sum 6^(j) and 6^(j) is 
free from e-dependent terms. 

The results of our numerical evaluation for the specific 
contributions for the 2P 1 / 2 state are 



M2P 1/2 ) = 0.33712, 
fa(2P 1/2 ) = 2.12732, 
0e(2Pi/ 2 ) = - 0.12830 . 
/3j(2P V2 )= - 0.52003 . 



and thus 



/3(2P 1/2 ) = 1.81611. 
For the 2P 3 / 2 state, we have 

f3 H (2P 3/2 ) = - 0.32557, 
(3^{2P m ) = 2.12732, 
(3e(2P 3/2 ) = - 0.12830, 
/3j(2P 3 / 2 ) = - 1.30008, 



and therefore 



/3(2P 3/2 ) = 0.37337. 



(68a) 
(68b) 
(68c) 
(68d) 

(69) 

(70a) 
(70b) 
(70c) 
(70d) 

(71) 



-(Za) 2 f3j(j), 

IT 



(65) 



where 5 jo is defined in Eq. (|23[) . Note that the correction 
to the current is ultraviolet finite and therefore does not 
contribute to the logarithmic part of the correction. The 



VII. VACUUM POLARIZATION 

The leading vacuum-polarization correction of order 
a(Za) 2 Ep is induced by a matrix element of the radia- 
tively corrected external magnetic field (corresponding 



10 



to a diagram with the hyperfinc interaction inserted into 
the vacuum-polarization loop). The corresponding cor- 
rection is (see, e.g., Ref. (25j ) 



t 2 



2t 2 



x / dr(l + 2mrt)e~ 2mrt g{r) f(r) 



(72) 



where g(r) and f(r) are the upper and the lower com- 
ponents of the reference-state wave function. To leading 
order in Z a, this expression is evaluated for P states to 
yield 




FIG. 2: (Color online) Same as Fig.Q] but for the 2P 3 / 2 state. 
The higher-order remainder function g 3 / 2 (Za) is plotted as a 
function of Z, with numerical values for g 3 / 2 (Za) given in Ta- 
ble H] The point at Z = is lim Za _o 53/2 (Za) = a 2 o(2P 3 / 2 )- 



tVP/„-\ - 01 ^2 2 " 2 - 1 



8C e t -{j) = -{Zaf~- 



(73) 



VIII. FINAL RESULTS 



Up to order a(Za) 2 Ep the vacuum-polarization correc- 
tion thus vanishes for P3/2 states. 



TABLE I: Higher-order remainder functions g 1 ^ 2 (Za) and 
gs/2 (Za) for the 2P l / 2 and 2P^/ 2 states, respectively, as ob- 
tained recently in Ref. [23[. The value of a employed in the 
calculation is a" 1 = 137.036, and the numerical uncertainty 
of the all-order (in Za) calculation due to the finite number 
of integration points is indicated in brackets. 



z 


g 1/2 (Za) 


g 3/2 (Za) 


1 


3.437410(86) 


0.12609(18) 


2 


3.369521(25) 


0.079405(55) 


3 


3.300616(14) 


0.032176(39) 


4 


3.231062(9) 


-0.015499(29) 


5 


3.161062(6) 


-0.063528(21) 


6 


3.090757(4) 


-0.111933(16) 


7 


3.020233(3) 


-0.160663(12) 


8 


2.949549(2) 


-0.209698(9) 


9 


2.878748(2) 


-0.259028(7) 


10 


2.807853(2) 


-0.308643(5) 




FIG. 1: (Color online) We show the higher-order remainder 
function g 1 / 2 (Za) as a function of Z. Numerical values for 
gi/ 2 (Za) are given in TableU The point at Z — is given by 
the coefficient a 2 o{2P 1 / 2 ). 



Summarizing our calculations, we conclude that the 
QED correction to the HFS of P states can be cast into 
the form 

AE nPj =E F [nP j )SUj) (74) 
= E F {nPj)- [aoo + {Za) 2 (a 2 i ln[(Za)" 2 ] + a 20 )] , 



which is valid up to order a(Za) 2 Ep. According to 
Eq. ([60|) . the results for the leading coefficients are 

a o(nPi/ 2 ) = 2 > a oo(nP 3 / 2 ) = -g • (75) 

The logarithmic part of the correction, as given by 
Eq. (021), is 



a 2 i(nP 



1 



1/2J 



o 31 (nP 3 / 2 ) =0. (76) 



In particular, there are no squared logarithmic terms. As 
follows from Eqs. (JBOJ), ||6T ]) . (|66 ]l . and fl75J), the total non- 
logarithmic contribution a 2 o is the sum of a self-energy 
(SE) and a vacuum-polarization (VP) correction, 



The results read 



a S 2 $(nP 1/2 ) 

a%o(nP 1/2 ) 

<4o(nP s/2 ) 
4^(^3/2) = 0. 



_ _ „SE 1 „VP 

a 2 Q — a 20 + a 20 



7 



(77a) 



19 

144 + 8^ + 16n 2 
+ P(nP 1/2 ), 

2 ti 2 -1 



In 2 



5 n 2 
109 



11 



1 



2880 32 n 30n 2 



(77b) 
(77c) 

■0(nP 3/a ), (77d) 
(77e) 
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Using the results for the (3 terms in Eqs. (|69|) and (|7Tj) . 
we obtain, in particular, 

4$(2P 1/2 ) = 3.50343 , (78a) 

c4 E (2P3 /2 ) = 0.17198. (78b) 

These results can be compared to numerical data at low 
Z presented in Ref. (23|. 

Indeed, for low-Z one-electron ions, a nonpcrturbative 
(in Za) calculation of the self-energy correction to the 
hyperfine splitting has recently been carried out [23[ . The 
numerical values for the self-energy correction AE„p. to 
the hfs find a natural representation as 

AE nP] = E F {nPj) 

x - [a 00 + (Za) 2 (a 21 ln[{Za)- 2 } + g 3 {Za))] , 

where gj(Za) is a remainder function which approaches 
the O2o coefficient for Za — > 0, 

lim gj(Za) = a,2o(nPj) . (79) 

Za— >0 

Numerical values for the remainder functions gj(Za) are 
given m Table [D In Figs. Q] and El we plot the higher- 
order remainder against the nuclear charge number Z. 
The numerical data are consistent with the next higher- 
order term in the expansion of AE n p j being a correction 
of order a(Za) 3 (no logarithm). 

IX. CONCLUSIONS 

The hyperfine structure of P states is an interesting 
physical problem. In accurate measurements of the clas- 
sical (2Pj-2S) Lamb shift in atomic hydrogen, both the 
hyperfine effects of P as well as of S states have to be 
carefully accounted for before a meaningful comparison 
of theory and experiment can be made. In alkali-metal 
atoms, the hyperfine structure of P states is also of great 
experimental interest [H, HH, H(| . 

In the present investigation, we analyze QED cor- 
rections to the hyperfine splitting of nP\j 2 and nP 3 / 2 
states in hydrogenlike systems, up to order a(Za) 2 Ep, 
where Ep is the Fermi splitting. Our calculation relies 
on the separation of the electronic from the nuclear de- 
grees of freedom as described in SecHH effectively reduc- 
ing the problem to an electronic self-energy type calcu- 
lation. The identification of the nonrclativistic degrees 



of freedom relevant to our investigation is accomplished 
by the Foldy-Wouthuysen transformation as described in 
Sec. IIIII A nuclear-spin dependent correction to the elec- 
tronic transition current is identified [see Eq. (|22p]. Wc 
show (see Sec. IIV|) that squared logarithmic corrections 
of relative order a(Za) 2 \n 2 [(Za)~ 2 ]Ep are completely 
absent for P states, whereas for P 3 / 2 , even the single 
logarithmic term of relative order a(Za) 2 \n[(Za)~ 2 ]Ep 
vanishes. This finding is interesting in view of different 
conjectures described in the literature [2lj . 

In order to address the nonlogarithmic correction of 
relative order a(Za) 2 Ep, we split the calculation into 
a high- and a low-energy part and match them via an 
intermediate overlapping parameter e that separates the 
scales of high-energy and low-energy photons (see Secs.lVl 
and lVip . This parameter is noncovariant but turns out to 
lead to a concise formulation of a problem which is oth- 
erwise rather involved. The high-energy part is treated 
in Sec. |V] and is seen to lead to form-factor type correc- 
tions. For the low-energy part treated in Sec. IVII the 
correction to the electron's transition current induced 
by the hyperfine interaction is crucial. This correction 
can be obtained via a Foldy-Wouthuysen transformation 
(Sec. IIIip . The discussion of vacuum-polarization correc- 
tions (see Sec. IVIip and a brief summary of the results 
obtained (Sec. IVIIip conclude our investigation. 

We reemphasize once more that the logarithmic a 2 \ 
coefficient vanishes for P 3 / 2 states [sec Eq. (|34p j. and the 
afo" self-energy coefficient [sec Eq. (|78bp ] also is numeri- 
cally small for 2P 3 / 2 . These two observations account for 
the numerically small results obtained in the all-order cal- 
culation [23| for the self-energy corrections to the HFS of 
this state. Indeed, the QED self-energy corrections to the 
hyperfine splitting of 2P 3 / 2 states are surprisingly small 
at low Z . This behaviour is naturally attributed to the 
less singular behaviour of the P 3 / 2 states at the origin in 
comparison to that of the Pi/ 2 states. 
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